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FOREWORD 

This  document  presents  the  results  of  work  performed 
by  personnel  of  the  Fluid  Mechanics  Section  of  Lockheed's 
Huntsville  Research  &  Engineering  Center.  This  document 
was  prepared  for  the  Engineering  Analysis  Division  of  the 
NASA- Johnson  Space  Center,  Houston,  Texas  under  Contract 
NAS9-13429,  Barney  B.  Roberts,  technical  monitor. 
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INTRODUCTION  AND  SUMMARY 


This  document  presents  a  description  of  the  modifications 
made  to  the  TRAN72  computer  program  to  meet  the  general  re¬ 
quirements  of  Lockheed’s  two-phase  flow  analysis  computer 
programs  (specifically  the  reacting  multi-phase  computer  pro¬ 
gram),  and  provides  instructions  for  operating  the  modified 
TRAN72  program. 


The  modified  TRAN72  program  is  available  for  external 
distribution. 
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Section  1 

CHEMICAL  EQUILIBRIUM  CALCULATION  MODIFICATIONS 

The  TRAN72  computer  program  (developed  by  NASA-Lewis  Research 
Center  (Ref.  1))  was  synthesized  by  combining  a  program  for  the  transport 
properties  calculation  with  the  CEC  71  program  (Ref.  2)  for  the  thermo¬ 
dynamic  properties  calculation.  The  TRAN72  program  was  subsequently 
modified  to  meet  the  requirements  of  Lockheed's  reacting  multi-phase 
computer  program,  RAMP  (Ref.  3).  The  requirements  satisfied  were: 

(1)  calculation  of  the  theoretical  rocket  performance  (for  both  equilibrium 
and  frozen  compositions)  during  a  "gaseous -only "  expansion,  after  a  two- 
phase  combustion  chamber  calculation;  (2)  inclusion  of  a  sudden  pressure 
freeze  during  a  two-phase  equilibrium  expansion  calculation;  and  (3)  auto¬ 
mated  communication  of  these  properties  to  the  RAMP  program. 

Modifications  were  made  to  the  TRAN72  chemical  equilibrium  calcu- 
lational  scheme  in  order  to  generate  thermochemical  data  consistent  with 
the  assumptions  utilized  in  the  RAMP  program  formulation.  The  assump¬ 
tions  being  addressed  in  the  RAMP  program  are; 

•  The  total  mass  of  the  mixture  is  constant. 

•  The  total  energy  of  the  mixture  is  constant. 

•  The  gas  obeys  the  perfect  gas  law  and  is  either  chemically 
frozen,  in  chemical  non equil ibrium  or  in  chemical  equilibrium, 

•  There  is  no  mass  exchange  between  the  phases. 

•  The  particles  are  inert. 

In  the  modified  TRAN72  calculational  scheme,  the  chamber  calculations 
are  performed  initially  with  the  condensed  species  considered.  The  total  mass 
and  total  enthalpy  of  the  mixture  are  then  adjusted  by  removing  the  mass  and 
enthalpy  associated  with  the  condensed  species  predicted  to  exist  in  the  chamber 
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alter  combustion.  The  total  mass  adjustment  is  made  by  removing  the  appro¬ 
priate  amount  of  mass  of  each  of  the  elements  which  comprise  the  condensed 
species  that  exist  in  the  chamber.  The  total  enthalpy  is  adjusted  by  removing 
the  enthalpy  associated  with  the  condensed  species  that  exist  in  the  chamber. 
The  adjusted  elemental  mass  balance  relationships  and  the  adjusted  total  en¬ 
thalpy  are  then  referenced  to  the  adjusted  total  mass  of  the  mixture.  All  con¬ 
densed  species  are  then  removed  from  the  list  of  possible  products  being 
considered  by  the  program.  The  chamber  calculations  and  subsequent  equi¬ 
librium  chemistry  expansion  are  then  made  with  a  gaseous-only  composition. 
Upon  completion  of  the  thermodynamic  calculations,  the  transport  properties 
are  calculated  in  the  manner  described  in  Ref.  1.  The  resultant  equilibrium 
chemistry  expansion  and  corresponding  transport  properties  data  are  for  the 
case  in  which  there  is  no  heat  transfer  between  the  condensed  and  gaseous 
species  during  the  equilibrium  chemistry  expansion  process.  To  account 
for  the  effects  of  the  heat  transfer  that  does  take  place  between  the  condensed 
and  gaseous  species  during  the  flowfield  calculations,  additional  thermo¬ 
chemical  data  is  required.  To  generate  the  required  data,  the  total  enthalpy 
of  the  gaseous -only  mixture  is  perturbed  (mass  is  held  constant)  and  the 
thermochemical  data  calculational  scheme  is  repeated.  The  total  enthalpy 
is  repeatedly  perturbed;  the  result  being  an  array  of  equilibrium  expansion 
processes  and  corresponding  transport  properties,  each  with  a  different 
degree  of  heat  transfer  between  the  two  phases. 

Fast  experience  in  thermodynamical  modeling  of  rocket  exhaust  flows 
has  indicated  that  many  chemical  systems  experience  a  transition  from  equi¬ 
librium  to  frozen  chemistry  during  the  expansion  process.  The  standard 
TRAN72  program  has  an  existing  option  to  treat  this  problem.  Under  the 
pressure  freeze  option  the  chamber  and  initial  expansion  calculations  are 
made  assuming  equilibrium  chemistry.  At  a  predetermined  pressure  ratio 
(chamber  to  local  static)  the  chemistry  of  the  system  is  frozen  and  the  re¬ 
mainder  of  the  expansion  is  completed  with  frozen  chemistry.  With  this 
option,  the  transport  properties  are  calculated  in  the  manner  outlined  in 
Ref.  1. 
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The  thermochemical  and  transport  data  are  communicated  to  the 
RAMP  computer  program  automatically  through  the  use  of  a  magnetic  tape 
(or  rapid  access  storage,  i.e.,  disc,  FASTRAN,  etc.).  Creation  of  the  data 
tape  (or  file)  is  accomplished  by  means  of  an  additional  subroutine  (MOCDAT) 
added  to  the  TRAN72  program.  Logic  is  provided  in  this  routine  for  creation 
of  a  new  data  tape  (or  file)  and  adding  data  to  an  existing  Master  data  list. 

Each  data  case  must  be  identified  with  a  unique  case  name  which  is  subse¬ 
quently  used  by  the  RAMP  program  to  determine  if  thermodynamic  data  are 
available.  An  additional  namelist  has  been  added  to  the  run  stream  to  control 
use  of  the  options  available  in  the  MOCDAT  subroutine. 

The  modified  TRAN72  program  is  used  to  generate  thermodynamic  and 
transport  properties  of  the  gaseous  phase  of  the  products  of  combustion  being 
considered  in  a  two-phase  flow  analysis.  Control  of  the  program  function  for 
this  application  is  handled  through  three  input  groups:  the  reactant  data  cards, 
the  $INPT2  namelist,  and  the  $RKT1NP  namelist.  A  detailed  description  of 
the  standard  TRAN72  program  input  is  given  in  Ref.  1.  The  thermodynamic 
data  required  for  this  application  are  calculated  using  the  RKT  option  under 
the  $INPT2  namelist.  Selection  of  this  option  permits  calculation  of  theoretical 
rocket  performance  for  both  equilibrium  and  frozen  compositions  during  ex¬ 
pansions.  The  variables  MOC2P,  PARTHT,  QDOTP  and  NQI  have  been  added 
to  the  $INPT2  namelist.  The  MOC2P  variable  controls  the  selection  of  the 
two-phase  flow  analysis  option  (MOC2P=T).  The  variables  PARTHT,  QDOTP 
and  NQI  control  the  selection  (PARTHT=T)  and  use  of  the  variable  total  en¬ 
thalpy  option  when  the  effects  of  heat  transfer  between  the  condensed  and 
gaseous  species  are  to  be  determined  in  a  two-phase  flow  analysis.  When 
PARTHT=T,  QDOTP  is  set  equal  to  the  amount  by  which  the  total  enthalpy 
of  the  gaseous -only  mixture  is  to  be  perturbed.  NQI  is  set  equal  to  the 
number  of  QDOTP  values  input.  The  specific  values  of  the  ratio  of  chamber 
to  local  static  pressures  (P  /P)  at  which  thermodynamic  and  transport  data 
are  generated  are  input  to  the  program  in  the  $RKTINP  namelist.  The  pres¬ 
sure  freeze  option  is  activated  by  setting  the  variable  NFZ  under  the  $RKTINP 
namelist  equal  to  the  number  of  the  pressure  ratio  at  which  transition  from 
equilibrium  to  frozen  chemistry  is  to  occur.  (The  chamber  is  considered  to 
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be  number  one,  the  throat  number  two,  etc,).  Freeze  pressures  may  be  the 
chamber  value  or  any  supersonic  pressure.  No  provision  is  made  for  freeze 
pressures  between  chamber  and  throat.  The  parameters  which  are  generally 
utilized  by  the  RAMP  program  are  local  Mach  number,  static  pressure  and 
temperature,  isentropic  coefficient  (gamma),  molecular  weight,  entropy, 
Prandtl  number,  viscosity,  specific  heat  at  constant  pressure  and  the  total 
enthalpy  (gas  only).  These  parameters,  with  the  exception  of  the  total  en¬ 
thalpy,  are  calculated  for  each  value  of  (P  /P)  ratio  by  the  program.  A 
detailed  description  of  the  logic  involved  in  the  standard  TRAN72  program 
computation  is  presented  in  flow  chart  form  in  Ref.  1.  This  information  can 
be  consulted  for  an  in-depth  understanding  of  the  calculational  scheme. 

To  automatically  create  a  tape  for  communication  with  the  RAMP  pro¬ 
gram  requires  that  one  of  the  two  tape-write  options  be  selected  (MOCT=T, 
or  MOCTF=T)  under  the  $INPT2  namelist.  The  MOCT  variable  is  utilized 
when  the  thermochem  data  are  to  be  run  completely  under  the  equilibrium 
assumption.  The  MOCTF  variable  is  utilized  when  the  thermochem  data  are 
to  be  run  completely  or  partially  frozen.  If  one  of  these  options  is  selected 
an  additional  namelist,  $TAPGEN,  must  be  input  to  control  the  tape-write 
function  and  the  input  of  the  case  name  card.  The  $TAPGEN  data  are  input 
after  the  $INPT2  data  but  prior  to  the  case  name  card  and  $RKTINP  namelist 
inputs.  Table  I  summarizes  the  program  variables  added  to  the  modified 
TRAN 7 2  program. 

Four  example  cases  showing  the  required  input  format  for  creation  of 
thermodynamic  data  for  typical  rocket  performance  problems  are  presented 
in  Table  2.  Case  1  is  the  required  input  to  perform  a  calculation  of  theoret¬ 
ical  rocket  performance  for  both  the  equilibrium  and  frozen  composition 
assumptions  during  an  isentropic  expansion.  (No  tape  is  generated.)  Case  2 
is  the  same  as  Case  1  except  that  a  tape  for  communication  with  other  pro¬ 
grams  is  generated  for  the  frozen  composition  assumption  during  expansion 
(MOCTF  =T).  Case  3  is  the  required  input  format  for  creation  of  thermo¬ 
dynamic  data  for  use  with  the  RAMP  program  (MOC2P=T);  a  tape  is  generated 
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Table  1 

ADDITIONAL  INPUT  VARIABLES  FOR  MODIFIED*  TRAN72  PROGRAM 


$INPT2  NAMELIST 

Variable 

Dimension 

Type 

Common 

Label 

Value 

Before 

Read 

Comment 

MOCT 

1 

L 

HREC 

F 

Selects  tape-write  option 
if  true  for  equil.  run  . 

MOC2P 

1 

L 

HREC 

F 

Selects  two-phase  flow 
analysis  option  if  true. 

MOCTF 

1 

L 

HREC 

F 

Selects  tape-write  option 
if  true  for  frozen  and 
pressure  freeze  options. 

PARTHT 

1 

L 

TWOPAS 

F 

Selects  variable  total 

enthalpy  option  if  true 
for  two-phase  analysis 

QDOTP** 

run. 

26 

R 

TWOPAS 

0.0 

Set  equal  to  the  amount 
by  which  the  total  enthalpy 
of  the  gaseous-only  mix¬ 
ture  is  to  be  perturbed. 

NQI 

1 

I 

TWOPAS 

0 

Set  equal  to  the  number 

of  QDOTP  values  input. 

$TAPGEN  NAMELIST 

IREAD 

1 

I 

1 

If  equal  0,  new  data  added 
to  master  data  tape  list; 
if  equal  1  data  written  on 
new  data  tape. 

IO 

1 

I 

— 

8 

Tape  unit  of  old  master 

1 

1 

*** 

IN 

1 

I 

- 

10 

tape  list. 

Tape  unit  of  new  data  tape. 

Case  Name  Card 

F  ormat: 

6A4 

Routines  modified  from  the  original  TRAN72  program  are:  LINK.MAINl 
REACT,  SEARCH,  EQLBRM,  ROCKET,  RKTOUT,  OUT1,  TRANSP,  OUT. 

The  values  of  QDOTP  must  always  be  input  in  ascending  order  (from  the 
most  negative  to  the  most  positive). 

*** 

When  running  multiple  cases,  the  data  of  the  last  case  must  always  be 
placed  on  tape  unit  10  if  it  is  to  be  communicated  automatically  to  the 
RAMP  or  VOFMOC  programs. 
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Table  2 

EXAMPLE  CASES  SllOWINC  THE  REQUIRED  INPUT  FORMAT 
FOR  CREATION  OF  THERMODYNAMIC  DATA  FOR  TYPICAL 
ROCKET  PERFORMANCE  PROBLEMS 

Case  1:  Required  input  to  perform  a  calculation  of  theroretical  rocket 
performance  for  both  the  equilibrium  and  frozen  composition 
assumptions  during  expansion.  (No  tape  is  generated.  ) 


mcc. r  an  r s 

M  QQ  1.00  o*o  r 

J.  t,oo  1.00  0*0 

(Insert  Blank  Card) 

NaMlL I  STS 
& 1 

KK  r  -  T  *  U-T  »nA^l  =  U0U01  «  •  00  «  oh  -  T  «i«ilX  =  o#C 

i»Li  *a 
i>Hrc  1  1  NP 

10.  .j0«  «  J  0  •  ,  i00.  .  O00  «  .  1000*  »  SO  00  .  .  J0000  »  «  IOOOoj.  . -JOOOOO. 

SL.aj 


Case  2;  Same  as  Case  1  except  that  a  tape  for  communication  with  other 
programs  is  generated  for  the  frozen  composition  assumption 
during  expansion  (MOCTF  =  T). 


KlAl T  mN I S 

M  c.uO  1.00  0.0  OCVb.lo  r 

oa.00  1.00  0*0 

(Insert  Blank  Card) 

N  IblS 

b  I  .  aHT  b 

kK  T  -  T  .  H  J  1  A  =  T  »  is.  A  Jt  -0000  1  «l-'  =  <i'OO.OO.Ol-  =  l  » >>.  iX-o.Ui  i'ijc  Tr  =  T 

it_'A  U 
b ThHoUv 

.  I  o  =  e  »  l  |M=  10 

it.AO 
CASL  z 
T 1 NH 

Kp-  1  o.  .  JO.  .-JO.  ,  too.  »J00.  .  1  000  *  vbOOO  •  .00000*  *  1  00000.  .  S00000  • 
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Table  2  (Continued) 

Case  3:  Required  input  for  <  reation  of  thermodynamic  data  for  use  with  the 

RAMP  program  (MOC2P=T);  a  tape  is  generated  for  t h<  equilibrium 
composition  assumption  during  expansion  (MOCT  =  T).  (The  effects 
of  heat  transfer  between  the  condensed  and  gaseous  species  are  not 
determined.) 


HEACTANTs 
AL.  1  .0 


C  6*004 

M 

1  0.009  J  .  <.  la 

N  *204 

FEd.O 

o 

3.0 

C  b*  lb 

H 

6.97  0  i . 1 7 

N  *0J 

N  1*0 

H 

A  . 0  J  4.0 

CL  l  *C 

(Insert  Blank  Card) 

OMIT 

AL  {  S  > 

AL  (L  ) 

OM  i  T 

ALNl S ) 

ALN 

OMI  T 

lCl3 

La-L4 

OM  l  T 

CH3 

CH4 

OMIT 

C2H6 

C302 

OM  J  T 

fl  {  S) 

Fl  (L  ) 

OM  l  T 

H20( 6 ) 

M20<  L ) 

NAME61STS 

’6  C 

1  •  OA 

.  4 

1  .96 
69.60 

Al.Ci_3  I  6  ) 

ALaCLb 

CH 

COU2 

CA 


0.0  629b. lb  F 

-12000.  629b. lb  F 

- 1 97300 .  6298.1b  F 
-28300.  S29b . 1 b  F 
-70690.  S29b . 1 b  F 


ACCL3<  L ) 

AL2Q2 

CM2 

C2CL2 

C5 

FECL2 ( a ) 


91NPT2 

HK  T  s  T  *P6  i  A  =  T  .  N  A  ;U  -OOuOi  .  P-bbA  •  00  .'■‘OCaP=  1  .t'ibCT  =  T 

*t.AO 
iT APOEN 

IPEAO-1 . Iu=8. (N-10 

*LNO 
CASE  3 
SkK T I NP 

PCP= 10. *30. .oo. . i JO . .600. .1000. .6000. • bCOOO . . 100000. • booooo. 
*END 


Case  4:  Same  as  Case  3  except  that  the  effects  of  heat  transfer  between 

the  condensed  and  gaseous  species  art  determined  (P\RTHT  =  T). 


REACTANTS 
AL  1.0 

16* 

0.0  S29U.15 

F 

C  6*064 

H 

10.089  o  .278 

N  •  264 

1<:*04 

-12000.  S298.1S 

F 

FE^.O 

0 

3.0 

•  4 

-197300.  S296.1S 

F 

C  6*  lb 

M 

6. *7  U  1.17 

N  .0  j 

1  •  Vb 

-2BJ00.  S296.15 

F 

N  1*0 

M 

* 

o 

o 

•r- 

c 

CL  l  •  0 

6V*e>0 

-70690.  S298.15 

F 

(Insert 

i'lan k  Card) 

OM  1  T 

At-(S) 

AL  (L  > 

ALl l  3 ( S ) 

ALCL3U-) 

OM  1  T 

ALN( 6 ) 

ALN 

AL2CL6 

AL202 

OM  I  T 

CCL3 

CCL4 

CH 

CM2 

I  T 

Ch3 

CH4 

C0CL2 

C2CL2 

OMI  T 

C2Hb 

C302 

C  4 

C5 

OM  1  T 

FE(SI 

FE  <  L ) 

FtCL2 l S  > 

FECL2 ( L. ) 

OM  f  T 

NAMEL I  STS 

H20C  S 1 

H20(Lt 

A  I NPT2 

WKT=1  *P->  lArT  ikftJt  '  OOOO  1  .  P  =  5bA.0‘J.  MOC2P-=  f  ".-iuCT  =T.PAPThT«T« 
ol>0TP  =  -600  .o.-300«0.0.0.100.0«NjI*a 

atNO 

1  I AP0tN 

i  we  a  j  s  i  ,io=a,i,\=io 

tlNb 

3.J0TP  Te^T  CAoL 

4WKT1NP 

PCP- 10.  <30.  •  >('..  100. .boo. -  I  000. .booo . .60000. . 1 orr  OO, , 500000* 
*tNO 
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for  the  equilibrium  composition  assumption  during  the  isentropic  expansion 
(MOCT-T).  (The  effects  of  heat  transfer  between  the  condensed  and  gaseous 
species  are  not  determined.)  Finally,  Case  4  is  the  same  as  Case  3  except 
that  the  effects  of  heat  transfer  between  the  condensed  and  gaseous  species 
are  determined  (PARTHT  =  T). 

2.2  USE  OF  THE  MODIFIED  TRAN72  PROGRAM  WITH  THE  STANDARD 
VOFMOC  PROGRAM 

The  TRAN72  program  has  been  modified  to  meet  the  requirements  of 
the  RAMP  computer  program.  The  data  tape  (or  file)  created  for  communi¬ 
cation  with  the  ramp  program  contains  additional  data  not  required  by  the 
VOFMOC  program  (Ref. 4).  For  that  reason,  the  tape  read  statement  and 
format  statement  in  subroutines  GASTAP  and  GASRD;  respectively,  must 
be  modified  to  read  the  additional  data  as  "dummy"  variables.  The  follow¬ 
ing  statements  must  be  changed  in  the  above  subroutines  before  the  data  tape 
generated  by  the  TRAN72  program  can  be  read  correctly  by  the  VOFMOC 
prog  ram. 

•  Subroutine  GASRD 

Old  Statement:  l  FORMAT  <  4A6 . 5X  ,A  3  .t>X  .  I  2  .  iiX  «  1  L  > 

New  Statement:  i  FORMAT  (  oA4  <  5X  *  A  i  *  6X  *  I  d. «  3X  «  1  £  ) 

•  Subroutine  GASTAP 

Old  Statement:  10  READ!  10)  <6ETA(  I  )  »  I  =  1 . 4  )  .  I  OF  .  I  i 
New  Statement:  10  READ  (  1  0  )  <  tJtTA  l  I  )  »  1  =  1 . 4  )  »OU  *  DU  *  I  OF  ♦  1  ^ 


There  are  no  other  limitations  placed  on  the  use  of  the  modified  TKAN72  pro 
gram.  Further  assistance  on  running  the  program  may  be  obtained  from  the 
authors. 
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CONCLUSIONS 

The  TRAN72  computer  code  (Ref.  1)  has  been  modified  by  Lockheed- 
Huntsville  personnel  to  provide  thermochemical  data  that  arc  compatible  with 
the  RAMP  computer  code  assumptions.  The  basic  capability  of  the  code  has 
been  left  undisturbed.  Thermochemi  1  composition  analyses  check  cases 
have  been  run  on  the  ABMDA.  7600  and  on  the  Univac  1108  Exec  8  at  Huntsville, 
Alabama. 

The  TRAN72  program  is  a  large-scale  scientific  endeavor  with  great 
capability  and  flexibility.  Programs  of  this  size  however,  are  rarely  finished 
but  usually  undergo  periodic  revision  to  increase  capability,  incorporate  sys¬ 
tem  changes  and  corrections  of  defects  which  might  occur  in  untested  situations. 
With  this  in  mind,  it  is  suggested  that  any  receiving  organization  address  ques¬ 
tions  concerning  the  basic  TRAN72  program  to  the  authors  of  Ref.  1.  Questions 
concerning  the  TRAN72  modifications  should  be  directed  to  the  Lockheed- 
Huntsville  authors  at  the  following  address: 

Lockheed -Mi  ssile  s  &  Space  Company,  Inc. 

Huntsville  Research  &  Engineering  Center 

P.  O.  Box  1103 

Huntsville,  Alabama  35807 
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